X-ray refractive index of laser-dressed atoms 
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We investigated the complex index of refraction in the x-ray regime of atoms in laser light. 
The laser (intensity up to lO^^W/cm^, wavelength 800 nm) modifies the atomic states but, by 
assumption, does not excite or ionize the atoms in their electronic ground state. Using quantum 
electrodynamics, we devise an ab initio theory to calculate the dynamic dipole polarizability and 
the phot oabsorpt ion cross section, which are subsequently used to determine the real and imaginary 
part, respectively, of the refractive index. The interaction with the laser is treated nonperturbatively; 
the x-ray interaction is described in terms of a one-photon process. We numerically solve the 
resolvents involved using a single-vector Lanczos algorithm. Finally, we formulate rate equations 
to copropagate a laser and an x-ray pulse through a gas cell. Our theory is applied to argon. 
We study the x-ray polarizability and absorption near the argon K edge over a large range of 
dressing-laser intensities. We find electromagnetically induced transparency (EIT) for x rays on the 
Ar Is 4p pre-edge resonance. We demonstrate that EIT in Ar allows one to imprint the shape 
of an ultrafast laser pulse on a broader x-ray pulse (duration 100 ps, photon energy 3.2 keV). Our 
work thus opens new opportunities for research with hard x-ray sources. 

PACS numbers: 32.30.Rj, 32.80.Fb, 32.80.Rm, 42.50.Hz 
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I. INTRODUCTION 

We study the interaction of atoms with two-color light. 
Specifically, we consider an optical laser with moderate 
intensity and x rays. In experiments this light would 
be obtained, for instance, from an amplified Tiisapphire 
laser system with up to lO^^W/cm^ at a wavelength 
of 800 nm. The x rays would be produced by a third- 
generation synchrotron radiation source like Argonne's 
Advanced Photon Source 

There are several ways for the chronology between the 
two light pulses. First, there is a pump-probe setting 
where the laser pulse precedes the x rays. Especially the 
combination of a weak laser as a pump and the x rays as 
a probe has received a lot of attention, e.g., Ref. 2. For 
higher laser intensities, ionization of atoms takes place, 
producing aligned, unoccupied atomic orbitals. They are 
probed by using xuv light or x rays to excite an inner- 
shell electron into them [1, 0, H, Q . Second, there is a 
pump-probe setting where the laser pulse succeeds the 
X rays. This situation has not received much of a focus 
(see, however, Ref. 0}. Third, there is a simultaneous ex- 
position of atoms and molecules to two colors. Recently, 
we studied molecules exposed to a laser with an inten- 
sity close to but still below the excitation and ionization 
threshold. If the molecule has an anisotropic polarizabil- 
ity tensor, they may be aligned along the linear laser 
polarization axis [i,ll,i,[i0, 11]; the X rays serve as an 
in situ probe of the rotational molecular dynamics. Go- 
ing a bit higher in laser intensity, without exciting and 
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ionizing, is possible for rare gas atoms. In krypton, laser 
dressing caused a pronounced deformation of t he p ho- 
toabsorption cross section around the K edge @, [ij- In 
neon we even found a strong suppression of the x-ray 
absorption on the Nc Is ^ 3p resonance [1, |Tl|, [isl IT^. 

The name electror nag netically induced trans- 
parency (EIT) [H, [M [13 for X rays was coined to 
describe the strong reduction of absorption in neon. It is 
related to EIT for optical wavelengths [ll[li, HO, HH HI 
by being essentially given by the same A-type three-level 
model. The only difference in the case of EIT for x rays 
is the linewidth of the second unoccupied level which 
is not approximately zero. Due to high laser intensity, 
transparency over a large number of wavelengths arises. 
The difference between the x-ray absorption spectra for 
laser-dressed neon and krypton is predominantly due 
to the ten times higher decay width of the K vacancies 
in krypton (Fis = 2.7 eV [23|, [23|) compared with neon 
(Fis — 0.27 eV [25]). K vacancies are produced by x-ray 
absorption and decay by x-ray fluorescence and Auger 
decay m. Another difference between EIT for optical 
wavelengths and EIT for x rays lies in the fact that 
refraction and dispersion by matter in the x-ray regime 
is generally tiny. 

EIT enables ultrafast control of the absorption and 
dispersion of a gaseous medium on a femtosecond time 
scale because the coupling laser, typically a Ti:sapphire 
laser system, can be combined with the sophisticated 
pulse shaping technologies available for optical wave- 
lengths [26|. These pulse shapes can be imprinted on 
X rays leading to ultrafast pulse shaping technology for 
X rays [1, [ll|, [l^, [13] • The ability to create almost ar- 
bitrarily shaped x-ray pulses opens up perspectives for 
the control of the dynamics of the inner-shell electrons 
in atoms. This way of shaping x-ray pulses is similar 
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to our recent study of x-ray absorption by laser-aligned 
molecules where x-ray pulses are shaped by controlling 
the molecular alignment [l^ . As the rotational dynamics 
of molecules takes place on a picosecond time scale, the 
x-ray pulse shaping is, however, done on a much slower 
time scale. Electron bunch manipulation techniques for 
short x-ray pulse generation are presently developed at 
synchrotron radiation facilities |27| . Our approach com- 
plements these efforts in a very cost-effective way. 

In this work, we would like to supplement our previous 
investigations with a detailed account of the refraction 
and the dispersion of x rays by laser-dressed atoms. Our 
theory is applied to argon. Also we demonstrate that 
there is appreciable EIT for x rays for argon. As there 
is no substantial EIT effect for krypton, this means that 
EIT for hard x rays is predicted for the first time (the 
K electrons of argon can be ionized with x rays of a 
wavelength of 3.9 A Such wavelengths are already 

useful to resolve coarse molecular structures via x-ray 
diffraction in molecular imaging; our work offers a way of 
temporal control of the x-ray pulses in such experiments. 

The paper is structured as follows. In SecHH we evolve 
our quantum electrodynamic formalism of Ref. [l^ to in- 
corporate rudimentary many-electron effects. We deter- 
mine the dynamic dipole polarizability for x rays and 
the x-ray absorption cross section in Sees. Ill Al III Bl and 
III CI on an ah initio level. Resolvents are evaluated in 
Sec. IIIDI with a single-vector Lanczos algorithm. Us- 
ing these quantities together with the classical Maxwell 
equations in Sec. Ill El we devise the complex index of 
refraction for x rays of laser-dressed atoms. The propa- 
gation of pulses through this new medium is treated in 
Sec. IIIFI Computational details are given in Sec. Illll for 
the results which are presented in Sec. IIVI Conclusions 
are drawn in Sec. |Vl Finally, secondary physical pro- 
cesses of the pulse propagation from Sec. Ill Fl are treated 
in the appendix. 

Our equations are formulated in atomic units [28| . 
The Bohr radius 1 bohr = 1 oq is the unit of length 
and 1 Iq represents the unit of time. The unit of en- 
ergy is 1 hartree = lE/j. Intensities are given in units 
of lE/i V^«o^ = 6.43641 x lO^^Wcm-^. Electric po- 
larizabilities are measured in 1 a§ E^^^ = 1.64878 x 
10^'*^ C^ TP? J^^. The dielectric constant of the vacuum 
is en = T-. 

II. THEORY 

A. Atoms in an electromagnetic field 

The time-independent Schrodinger equation of the 
field-free atom is 

i?AT ^'f (n CTi, . . . , az-, t) = Es *f (n ai, . . . ,fz crz -.t) . 

Here, -STat is the electronic Hamiltonian of the atom 
which contains the Coulomb interaction of the Z atomic 



electrons with the nucleus and the two-particle inter- 
action among the electrons psi . [29| . The wave func- 
tions ^'f (ri CTi, . . . , rz CTz, with eigenenergy Es repre- 
sent the ground state (s — 0) and core-excited excited 
states (s > 1) of the atom [301 for the uth electron at 
position r\j with spin projection quantum number Uu 
for u g {1, . . . , Z}. 

The free electromagnetic field is represented by the 
Hamiltonian 
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k ^=1 

H„.,t^ = ^-ai a- -TV- cj- , (2b) 

EM,fe,A k k,\ k,\ k,X k ' ^ ^ 

where at and a- denote creation and annihilation op- 

k,X k.X 

erators, respectively, of photons in the mode k, A with 
wave vector fc, polarization index A, and energy [3l| . 
Here, indicates the number of photons in the 

mode fc, A of some initial state. The energy of the ini- 
tial state, Eq. (fTU)) in our case, has been adjusted such 
that the field energy is zero. This simplifies our equa- 
tions notably. However, physical quantities are of course 
independent of the absolute energy shift in Eq. 

We consider in this paper only two modes for the elec- 
tromagnetic field denoted by L and X for the laser and 
the X rays, respectively. The eigenstates of the free elec- 
tromagnetic field Hem are Fock number states 

|iVL)®|7Vx)-(a[)^M4)^^ |0) - (3) 

with the vacuum state 1 ) . 

The interaction of electrons with light is described by 

using the principle of minimal coupling to the electro- 
magnetic field [31I, [121 • The field is represented by the 
vector potential A{r) for which we assume the Coulomb 
gauge. The electrons are created and annihilated by the 
field operators ipl-{f) and ipai'r), respectively f29j. We 
assume that both laser and x-ray wavelengths are suffi- 
ciently large for the electric dipole approximation to be 
adequate. Then, the vector potential becomes indepen- 
dent of f. 

In dipole approximation, the form ^ of the interaction 
Hamiltonian with the x rays reads 

^ il j + ^^x] Mr) , (5) 

using the following mode expansion for the quantized vec- 
tor potential of the x rays |3l| : 



^x = \/-f7 [ex ax + ex a]^] . (6) 
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It is frequently referred to as velocity form [33| . 

In laser physics, one typically transforms the interac- 
tion 111 in dipole approximation to the so-called length 
form. This can, of course, be accomplished either for the 
Hamiltonian with semiclassical electromagnetic fields ^] 
or, as in our case, in the quantum electrodynamic frame- 
work [3l] |. Then, the interaction Hamiltonian with the 
laser field becomes 

Hi,L= Ji'U^r-ELMf^d^r, (7) 

where we use the quantized mode expansion 



Et 



(8) 



for the electric field of the laser mode [31| . The influence 
of the two-color light ^ in dipole approximation has 
been decomposed into Hi = Hil + Hix [Eqs. (O and 
©]. 

The full Hamiltonian of an atom in two-color light H = 
Hq + Hi consists of a strongly interacting part [12l | 



Hq = Hq i^ + Hq^x , 
Hq^l — Hat + Hem,l 
Hq^x = Hem,x , 



(9a) 
(9b) 
(9c) 



and a weak perturbation due to the interaction with the 
xrays. Hi = Hi,x M- 



B. Dynamic Stark effect 

We would like to determine the energy of the if-shell 
electrons in the two-color field. It is given by a pertur- 
bative expansion with respect to the x-ray field using 
a basis of laser-dressed energy levels of the atom. The 
initial state of the atom in light is given by the direct 
product of the electronic ground state wave function and 
the Fock states for the laser and the x-ray mode 



l*f>®I^L) 

\i)(g)\Nx) . 



(10a) 
(10b) 



The energy of the initial state is £"7.0 = (-^1 ^^o 1-^)- We 
assume that the electronic ground state of the atom is 
not noticeably altered by the light. Core-excited states 
with a vacancy in the K shell serve as final states for 
x-ray absorption and emission. We define for s =/= 0: 



I 



)|7Vx±i; 



(11a) 
(lib) 



Let us solve the strongly interacting laser-only prob- 
lem (I9bp first in terms of a non-Hermitia n, comp lex- 
symmetric representation of the Hamiltonian [sj, [s^ , 
using the laser-only basis (jllap for s, s' ^ 0: 



' s,i.i:s' .fj,' 



{^s,^,\Ho.L\^s'.,^^' 



(12) 



The eigenstates \F) of Eq. p2|) form the laser-dressed 
atomic energy levels with complex eigenenergies Ep. 

The interaction with the x-ray field is treated as a per- 
turbation of the laser-free atomic ground state; the im- 
pact of the laser dressing arises exclusively in the final 
states. The full basis states that include the x rays are 
defined in Eq. (|llb[) : the laser-dressed atomic levels are 
given by the direct product states | F"^ ) = | -F)'?' | NxH) 
which correspond to the energies Ep — ( F^ \ Hq \ F^ ) = 
Ep ± ujx- We use non-Hermitian perturbation the- 
ory [12, 37] to determine the complex energy 



Ei-cs — Ej^ — i r/2 



(13) 



of the initial state (|10b[) , the so-called Siegert energy [3J, 
,38] . The real part is the energy shift of the atomic level 
due to the laser and x-ray fields [13, HI] whereas T stands 
for the transition rate from the initial state to Rydberg 
states or the continuum. The energy is obtained 
from 




I\Hi\F^){F'\Hi\I) 



(14) 



E 



1,0 



El 



P2(wx) 

Assuming A'x S> 1, the interaction matrix ele- 
ments {I\ Hi\F^) and {I\Hi\F^) arc approximately 
the same. 

The first order correction Pi(cjx) in Eq. (|14p is found 
by combining Eqs. (HI), and ((TI 



Piiiux) = - {I\ 
ttZ 



E 



Vujx 



i2Nx + l) 



27r 



Vujx 



NxZ 



(15) 

The one-particle operator hereby acts on all Z elec- 
trons in the initial state (jlObp . The tiny correction y-^ 
after the second equals sign vanishes upon taking the 
limit \^ cx) in the end. The real quantity Pi{ujx) in 
Eq. (fT5|l is the energy shift due to the scattering of x rays 
and represents the ac Stark shift [37, 40]. We find the 
corresponding dynamic polarizability [40] from 



Re Pi{LOx)^-^ai{cJx) Flo 



(16) 



The peak electric field of the x rays is Ex,o = 
^/Snalx [33I . I4H for the x-ray intensity Ix- From 
Eqs. ([TS]) and (fT6|) . we obtain the first order atomic po- 
larizability 



ai(wx) 



Z 
7T2" 



(17) 
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In analogy to Eq. ([16]), we define the second order con- 
tribution to the dynamic polarizabihty [4fl|]. With the 
real part of P2{uJx) in Eq- p4|) . the well-known Kramers- 
Heisenberg form 

4 

a2(wx) = --p^ ReP2(wx) 



x.o 



ReJ2{I\H,\F+) {F+\H^\I) 



x,o 

^ 2 {Ep — Ej^q) 
{Ep — EifiY ~ 

(18) 

results ^,\^. 

Using the first and second order atomic polarizabilities 
from Eqs. p7p and (|18p . we can define the total atomic 
polarizabihty 



ai{ujx) = ai{ujx) + a2{ujx) ■ 



(19) 



We derive explicit equations for the polarizabihty a2{uJx) 
in the next section. 



C. Hartree-Fock-Slater approximation 

We solve the Schrodinger equation ^ with the 
help of the independent-electron approximation (l^ . [28l . 
[29} . Then, the Z-electron ground-state wave function 
is given by a Slater determinant of one-electron or- 
bitals <1>Q (r*! (Ti, . . . , fz crz)- Using the formalism of sec- 
ond quantization, we have 



1$^ 



n 



n,l,m,^ n,l,m,\. 



|0) (20) 



n^l,m occupied 



for a closed-shell atom. Here, ^„ ; ^ | and ; „ are 
creators of electrons in the spatial orbital (^„,;.,„(r) with 
spin projection quantum number up (ti = ^ = t) and 
down (cr = — i = I), respectively. The spatial orbitals 
are eigenfunctions of the one-electron Hamiltonian 



Subsuming these contributions, we obtain the effective 
one-electron atomic Hamiltonian [12] which we use in- 
stead of -ffAT in Eq. ((9b)) . 

With the help of the spin orbitals, we can expand the 
field operators 29] in the equations of Sees . Ill Al and [HBl 
as follows 



l,m,(T ! 



(22a) 



^l{r) = (22b) 

n,l,m 

Based on the independent-particle ground state (1201), w e 
construct the spin-singlet i^-shell-excited states [2^, [23] : 



1 

V2 



|^n,Z,m,/x) — /77 [^i,i,m,f ^1,0,0, T 



® |iVL-/i> 



1$ 



Nx±l] 



(23a) 



(23b) 



as basis states [see Eq. (fTTj) ]. Note that to is a conserved 
quantum number because the laser is linearly polarized; 
thus the matrix (fT2|) blocks with respect to to. Following 
Ref. m, the Kramers-Heisenberg form (|18p becomes 



a2(wx,'i? 



LX J 



Re 



1 

E 

m— — 1 



><m(^?Lx) 



p. — Els — ^X 



(24) 



Here, we omit the emission term "+" which results from 
decomposing the denominator of Eq. fTS]) to emphasize 
the similarity to the expression for the cross section. The 
initial-state energy is the iC-shell energy EpQ = Eis- 
The angular dependence of the polarizabihty (|24p is given 
by >fm(i9Lx) which is equal to cos^ t^lx for to = and 
equal to i sin^ ??lx for to = ±1. The radial dipole matrix 

elements are given by [l^ . 



HFS 



Uhfs(»') 



(21) 



They are characterized by the principal quantum num- 
ber n, the angular momentum ^, and the magnetic quan- 
tum number m (29j . To determine the effective one- 
electron central potential Vhfs ('') ) we made the Hartree- 
Fock-Slater mean-field approximation [4^. |43|. 

When X rays are absorbed, an electron may b e ej ected 
into the continuum leaving a core hole behind [4j. We 
use a complex absorbing potential (CAP) to handle such 
continuum electrons. The CAP is a one-particle operator 
which is added to ifups- It is derived from smooth ex- 
terior complex scaling [l^, IH, IH, \^ . Additionally, one 
needs to allow for the relaxation of core holes by x-ray 
fiuorescence and Auger decay P, HI] with a decay width 
of Pis- According to Eq. (fT5|) . this is accounted for b:v 
adding — iris/2 to all energies of core-excited states [1^ 



D. Lanczos solution of the polarizability 

The full diagonalization of the interaction of the atom 
with the laser to obtain the new energy levels of a laser- 
dressed atom becomes impractical for larger numbers of 
atomic orbitals and photons. To be able to solve the 
equations for the polarizability (|24p and the cross sec- 
tion [13] with acceptable computational effort, we em- 
ploy a single- vector Lanczos algorithm [i^ [s^]- A real 
symmetric version has been implemented by Meyer and 
Pal [m ; it was extended to complex symmetric matrices 
by Sommerfeld et al. [s^ . 

The expression for P2{ijJx) in Eq. (fTi)) can be rewritten 



P2(wx) = {I\HiP 



1 



Els — Ha 



PHi\I) 



(25) 
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with the projector on the core-excited basis states ()23bp : 



P- = 



E 



Let us define the resolvents 



(m),± 



n.l.m 



1 



,m,/i I 



(26) 



n.ljni.i^L I 



Els — Ho 



I u 



(27) 

the two cases for x-ray absorption and emission are only 
distinguished by a ^uiy^ term in the denominator. Apart 
from this, the expressions involves only the interaction 
with the laser. In fact, the resolvents P7)) can be ap- 
proximated in terms of a single-vector Lanczos run for 
the matrix [Eq. (fT2)) ]. To show this, let the com- 

ponents of the start vector [s^, [Mj be 



(m) 



I) 



(28) 



In Floquet approximation, the results are the same for 
the two cases and "— " of x-ray emission and absorp- 
tion; we use the former case throughout. Using Eqs. (pS)) . 
(P7)) . and (pS)) . expression (^5]) . is recast into 



-1 



(29) 



index of refraction. From the Maxwell equations, we ob- 
tain the equation for the electric field (ssl . [53| of a wave 
propagating in the z direction: 



d^E{z,t) 1 d^E{z,t) _ Att d^P(z,t) 



dz^ 



a<2 



(31) 



Let E{z,t) = E+(z,t) + E-{z,t) and E-{z,t) = 
[E+{z,t)]* [3l,|4l| with 



E+{z,t) = ii^oe'C^^-'^^*) . 



(32) 



Here, k is the wave number in the medium. The complex 
polarization P{z,t) describes the linear response of an 
atom to an electric field in terms of the complex electric 
susceptibility xi^x) via 



P±(z,t) ^ l.x{LOx)E^iz,t) . 



(33) 



We define P{z,t) = P+(z,t) + P"(z,<), P^{z,t) = 

[P+(z,i)]*,andPo = 47X(wx)i?o. 

Inserting P{z,t) and E{z,t) into Eq. ([?!]) yields the 
relation between the refractive index and the polariza- 
tion [H: 



"^(wx) = 1 + xi^x) 



2^2 



k'^C 



info 
Eo 



(34) 



After iVitor Lanczos iterations, we have reduced H, 



(m) 
0,L 



to the tridiagonal matrix T^") = Q(")Tjjr(m)g(m) 

where Q'™-* is the matrix of Lanczos vectors. The ma- 
trix t'"*-* has the eigenvalues A^™-*, and the eigenvec- 
tors i.e., we have x('")T2.(m)^(m) ^ ^(m)^ 

insert (g('")x('"))(Q(")x("))T = l on the left hand 
side of t/'™-* in Eq. (p9|) and its transpose on the right 
hand side of i/'™''^. The start vector ([28| is normal- 
ized at the beginning of the Lanczos iterations. Hence 

that g(™)T^m) ^ |j^(m)|| ^i^lj ^ (1,0,..., 0)T. 

Finally, Eq. (^5)) reads in the Lanczos basis 



P2(C^X) 



Niter 

E 

p=i 



A^ 



(m) 
PP 



-fcls - t^X - App 



(30) 

From P2('-^x) in Eq. (j30p . we obtain the dynamic atomic 
polarizability fTS]) by taking the real part of P2(a;x). The 
photoabsorption cross section follows from the imaginary 
part of P2(wx) using the prefactors in Ref. IT^ . 



E. Index of refraction 

We characterize the impact of a gaseous medium of 
laser-dressed atoms on x rays in terms of the complex 



Consequently, the propagation through the medium is 
described by Eq. PT|) in terms of 



E+{z,t) = i gi ["(c^x) fco^-<^xt] 



(35) 



with the wave number in vacuum /cq = ^ . Equation (|35p 
is used to form a solution E{z, t) — E^{z, t) + [E^{z, t)]* 
of the wave equation pT|) . i.e., Eq. ([34]) is satisfied upon 
letting k — n{u!x) ko- 

To make the connection of our classical electrodynamic 
equations to our quantum electrodynamic results from 
Sec. IIIB| we use the dynamic polarizability ([T9|) to ex- 
press the polarization (|33p of the medium with number 
density due to the x rays, 



Po = n# ai{ujx)Eo 



47r fco 



Eo , 



(36) 



where we add an imaginary intensity absorption term. 
It contains the absorption coefhcient fi = cr(ujx_) 
which involves the x-ray absorption cross section of the 
atom a{iux) Using Eqs. ^ and the index of 
refraction becomes 



n(wx) w 1 + 



Xi^x) 



l + 27rn#ai(u;x)+i;r7- , (37) 
Z ko 



where the square root of the right hand side of Eq. ([341) 
was approximated. The imaginary part in Eqs. p6p and 
(157)1 leads together with Eg. ([55)1 to an exponential decay 
following Beer's law (4ll. l4a|. 
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Photon energy [eV] 



F. Pulse propagation 

We assume copropagating laser and x-ray pulses which 
pass through a gas cell of length L that is filled with a gas 
of pressure p and temperature T (see Fig.|9]). To describe 
the evolution of the pulses, we employ the rate-equation 
approximation [4l| . The rate equations are solved on a 
numerical grid of A^grd points to discretize the distance 
covered during the propagation time interval. Let the 
laser radiation of angular frequency be linearly polar- 
ized along the z axis. The laser intensity is given by /l {t) . 
It is assumed to be constant perpendicular to the beam 
axis. We assume an x-ray pulse with a Gaussian enve- 
lope with peak flux Jx,o and a full width at half maxi- 
mum (FWHM) duration of rx: 

JxW = Jx.oe""'"'^*^' . (38) 

The peak flux of the pulse (at t = Q) depends on the 
number of photons per x-ray bunch like 

Jx.o = 2./i^^X(0). (39) 
V TT rx 

The factor X(0) — represents the peak of a Gaus- 
sian radial profile of a FWHM width of px- 

We assume that the laser light and the x rays both 
propagate with the speed of light. In the appendix, we 
made a number of estimates of the impact of secondary 
physical effects in the gas cell on the phase and group 
velocities of the laser and the x rays. 

We describe the interaction of the two-color light with 
the gas by the following rate equation for the x-ray pulse 

AJx(^j+i) - -n'^{zj)<y'{z,) J^{z,) Az , (40) 

which is evaluated for the individual steps j G 
{!,..., A^grd — 1} to move the pulse over the grid. Let 
a'{zj) = a{Ii,{zj / c)) and Jyi{zj) = JyJyZjjc). Further, 
n'^iyZj) represents the number density of argon atoms 
at Zj . Initially, it is n# inside the gas cell and zero other- 
wise. Equation (|40)) can be solved analytically in the gas 
cell; one obtains an exponential decay following Beer's 
law [4l|, for each infinitesimal section of constant in- 
tensity of the x-ray pulse. 



FIG. 2: (Color online) X-ray absorption cross section of argon 
for parallel laser and x-ray polarization vectors. In the upper 
panel, the intensity of the dressing laser is lO^'^ W/cm^; it is 
10^^ W/cm^ in the lower panel. 



III. COMPUTATIONAL DETAILS 

The computations presented in the ensuing Sec. IIVI 
were carried out using the dreyd and the pulseprop 
programs of the fella package 1541. We developed 
DREYD and PULSEPROP for Refs. [gI. I 111. [Tsl H For this 
paper, we added the treatment of the velocity form ^ of 
the interaction Hamiltonian with the x rays to dreyd. 

The computational parameters of dreyd are speci- 
fied following Ref. 12,. To solve the atomic electronic 
structure problem, we use the Hartree-Fock-Slater code 
of Herman and Skillman [5^ setting the Xa parame- 
ter to unity. The radial part of the atomic orbitals 
is represented on a grid of a radius of 60 ao using 
3001 finite-element functions, considering angular mo- 
menta up to Z = 7. We represent the radial Schrodinger 
equation in this basis set. From its eigenfunctions, we 
choose, for each I, the 100 lowest in energy to form atomic 
orbitals [l^ . The smooth exterior complex scaling com- 
plex absorbing potential is formed using the complex 
scaling angle 9 — 0.13rad, a smoothness of the path 
of A — 5aQ^, and a start distance of tq = 7 uq. We 
think of a Ti:sapphire laser system as a potential dress- 
ing laser which emits light at a wavelength of 800 nm, 
i.e., a photon energy of cjl = 1.55 eV. To converge the 
x-ray absorption spectra, we accounted for the emission 
and absorption of up to 12 laser photons in the Floquet- 
type matrix (H^). We set the K edge of argon to its 
experimental value of Eis = 3205.9 eV [ij as well as the 
linewidth of a, K vacancy Fis = 0.66 eV [HB]- Finally, 
we need to carry out A'itor = 4000 Lanczos iterations to 
converge P2{ujx) in Eq. ([50)1 . 

The PULSEPROP code assumes a gas cell with a length 
of 6 mm (Fig. [9]). It is filled with argon gas of a pres- 
sure of p = 5 atm at a temperature of T = 300 K. The 
number density of argon atoms follows from the ideal gas 
law rt^ = ^ = 1.22 X lO'^" cm~^. The x rays are tuned to 
the Ar Is ^ 4p resonance which has in our computations 
an energy of 3203.42 eV (peak position in Fig. [T|). Each 
x-ray pulse (|38p comprises nph = 10^ photons and has 
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FIG. 3: (Color online) X-ray absorption cross section of argon 
for perpendicular laser and x-ray polarization vectors. The 
dressing laser has an intensity of lO^^W/cm^ in the upper 
panel and of 10^^ W/cm^ in the lower panel. 
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a FWHM duration of rx = 100 ps. The x-ray beam is 
circular with a FWHM focal width of 10~^ cm. This im- 
plies a peak x-ray flux p9|l of Jx,o = 8.3 x 10^ ps~^cm~^. 
The laser pulse envelope has the shape 

4 

Il{zo,zl,z) = /l,o ^ cos(^^Tnj 

m=0 (41) 
^ g-41n2 [(2:-Zo+2(m-2)2:L)/zL]=^ _ 

The FWHM length of the constituting five individual 
laser pulses is defined by zl = ctl for the FWHM pulse 
duration tl = 1 ps. The pulse train is centered at zq. 



IV. RESULTS AND DISCUSSION 

In Fig. [U we show the x-ray absorption cross section 
of argon near the K edge without laser dressing. Clearly, 
the isolated Ar Is ^ Ap pre-edge resonance is discernible 
at 3203.42 eV. The much weaker Arls 5p pre-edge 
resonance can also be seen at 3204.72 eV. The contin- 
uum stretches beyond the K edge of argon which lies 
at 3205.9 eV ^. 

We expose the argon atoms to a linearly polarized 
800 nm laser with intensities up to lO^^W/cm^. The 
intensity remains well below the appearance intensity 
of Ar+ [H^l- This is an important assertion because our 
theoretical description does not account for ionization 
without prior x-ray absorption. 

In Fig. [21 we display the x-ray absorption cross section 
with laser-dressing and parallel laser and x-ray polariza- 
tion vectors. The cross section is shown for the two laser 
intensities lO^^W/cm^ and lO^'^W/cm^. The spectrum 
for lO^^W/cm^ in the upper panel of Fig. [5] resembles 
the one without laser in Fig. [TJ Yet the Ar Is Ap reso- 
nance is split into two close- lying peaks. Also small new 
bumps emerge in the spectrum, e.g., around 3206 eV. For 
lO^^W/cm^, the trends become more pronounced; the 
absorption spectrum is significantly modified compared 



FIG. 4: (Color online) Dependence on the intensity of the 
dressing laser of the x-ray absorption cross section of argon on 
the Ar Is 4p pre-edge resonance at 3203.42 eV for parallel 
laser and x-ray polarization vectors (solid black curve) and for 
perpendicular vectors (dashed red curve). The dash-dotted 
green line represents the cross section without laser dressing. 



to the laser-free case. The Ar Is ^ 4p resonance is split 
into two peaks which are separated by almost 2 eV. 

When one has tuned the x-ray energy to the peak of 
the Ar Is Ap resonance without laser dressing, then, 
turning on the laser, leads to a substantial suppression 
of x-ray absorption. This phenomenon was analyzed by 
us in detail in Ref. jl^ We termed it electromagnetically 
induced transparency (EIT) for x rays. The mechanism 
which leads to EIT for x rays was analyzed in Ref. [l^ in 
terms of a quantum-optical three-level model fl^. The 
model reproduces well our ab initio data proving that 
our interpretation of the phenomenon is adequate. The 
interpretation for argon is very similar. The model is 
formed by the levels Ar Is, Ar ls~^4p, and Ar ls~^4s. 
The X rays couple Ar Is and Ar ls~^4p whereas the laser 
couples Ar ls~"'^4p and Ar ls~^4s. The mechanism of EIT 
for X rays is very similar to EIT in the optical domain, the 
only difference being an appreciable decay width of the 
third state (Arls~^4s). The transparency is therefore 
due to a splitting into an Autlcr-Townes doublet and not 
primarily due to destructive interference of two excitation 
pathways ^21|] . 

In Fig. [31 we show the impact of laser dressing for 
perpendicular laser and x-ray polarization vectors. In 
contrast to the case of parallel polarization vectors in 
Fig. [21 the suppression of x-ray absorption is apprecia- 
bly smaller. In fact, the curve for the lower intensity 
in the upper panel looks nearly identical to Fig. [1] The 
Ar Is ~> 4p resonance in the lower panel is split into 
two peaks. Also minute bumps can be perceived beyond 
the K edge. Let us analyze the underlying quantum 
optical mechanisms: we have no direct dipole coupling 
of Arls~^4p and Arls~^4s states. This excludes the 
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FIG. 5: (Color online) X-ray polarizability of argon near the 
K edge without laser dressing. 



mechanism discussed for the case of parallel polariza- 
tion vectors. The suppression of the absorption there- 
fore is attributed to arise partly through a line broaden- 
ing caused by the laser dressing Additionally, the 
Ar ls~^4p and Ar ls~^3(i states are in resonance within 
the linewidth broadening due to the inner-shell decay. 
The radial dipole coupling matrix clement between the 
two states is even larger than the coupling between the 
Ar ls^^4p and Ar ls~^4s states. We conclude that these 
two states together with the ground state Ar Is form a 
three-level ladder-type model. Such models also exhibit 
EIT — although not in a strict sense (see Ref. [2l| and ref- 
erences therein for details) — which we observe here for 
the first time in the x-ray regime. 

In Fig. m we show the dependence of the x-ray ab- 
sorption cross section of argon on the Ar Is — > Ap reso- 
nance on the intensity of the dressing laser. The cross 
section without dressing laser, 128.6 kb, is indicated. 
Above a laser intensity of 10^^ W/cm^, the cross section 
drops steeply towards the value 35.6 kb at lO^^W/cm^ 
for parallel laser and x-ray polarization vectors. This 
corresponds to an overall reduction of the cross sec- 
tion with laser and without of a factor of 3.6. For 
perpendicular polarization vectors, only a reduction by 
a factor of 1.9 is found. We observe that the speed 
of the drop of the cross section decreases noticeably 
near IQ^^W /cm^. Eventually, the cross section seems to 
saturate beyond 10^'^W/cm^ before appreciable strong- 
field ionization of the ground-state atoms takes place 
around 10^'*W/cm^. 

These observations can be understood in terms of the 
Rabi frequency f^L — di?L,o of the laser coupling of the 
Ar4s and Ar4p levels. Here, d = {As\z\Ap) = 3.5 oq 
is the dipole coupling matrix element between the two 
states and E'l.o = \/8tt a Ix,q the peak electric field. 
For the laser intensities lO^^W/cm^, lO^^W/cm^, and 
10^^ W/cm^, we find the Rabi frequencies WE' ^^'^ 
2Q[s - The frequency fl^ should be comparable with the 
inverse core-hole lifetime of argon of 1 fs [s^ such that 
it beats the decay and absorption is suppressed. The 
r^L also determines the magnitude of the splitting of the 
lines in the Autler-Townes doublet [l3|- However, the in- 
creasing laser intensity also leads to a broadening of the 
Rydberg states due to ionization. The line broadening 
competes with the Rabi flopping and finally slows down 
the drop of the cross section with increasing laser inten- 
sity in Fig. m It potentially leads to a halt of the drop 




X-ray photon energy m,^ [keV] 



FIG. 6: (Color online) Refractive index of argon without laser 
dressing. The theoretical data are represented by the dashed, 
red curve. Semiexperimental data (using the Kramers-Kronig 
relation to determine the real part of the refractive index 
from the experimental cross section data) was taken from 
Fig. 4 in Ref. Is^ . The digitized experimental curve was fit- 
ted with n(a;x) - 1 = ao/^x' yielding ao = -672233 eV"i 
and ai = 2.096 for the two parameters. Then, the curve was 
shifted such that n(606 Ry) - 1 = -0.503 x 10* (see "This 
paper" column of Table V in Ref. [s^ . 



before ionization kicks in. 

The real part of the energy shift in the x-ray field de- 
termines the atomic polarizability (fT9|) . It is plotted for 
argon around the K edge in Fig. [51 Clearly the dis- 
persion curve of the Ar Is —> Ap pre-edge resonance is 
perceivable. 

From the atomic polarizability, we determine the re- 
fractive index at normal temperature and pressure [60l | 
from Eq. ([57)) . In Fig. [51 we compare our theoreti- 
cal results with the data from Fig. 4 of Liggett and 
Levinger [sl] . We find a very good agreement above the 
K edge of argon. Below 300 eV, the refractive index is not 
even qualitatively reproduced because our approximation 
breaks down that the refractive index is determined only 
by K-sheW electrons; the electrons in higher lying shells 
become significant. 

The atomic polarizability and thus the refractive in- 
dex is changed by laser dressing. In Fig. [71 we investi- 
gate the impact of parallel laser and x-ray polarization 
vectors and in Fig. [H we plot the dynamic polarizabil- 
ity for perpendicular polarization vectors. Generally, the 
polarization — and thus the refractive index ([37]) — is dom- 
inated by the A\ term. Hence the impact of the laser 
dressing is expected to be small because it is mediated 
by the p ■ Ay^ term. This expectation is confirmed by 
Figs. [3 and [51 the magnitude of refraction and dispersion 
of X rays is small with or without laser dressing. The 
shape of the curves in the two figures is, however, very 
different. While for parallel polarization vectors a sub- 
stantial change of the shapes is obtained for 10^^ W/cm^, 
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FIG. 7: (Color online) X-ray polarizability of argon for 
parallel laser and x-ray polarization vectors. In the upper 
panel, the intensity of the dressing laser is 10^^ W/cm^; it is 
10^^ W/cm^ in the lower panel. 



the corresponding curve for perpendicular polarizations 
is nearly unchanged compared with Fig. [5l The different 
behavior again reveals the different underlying quantum 
optical process already discussed for absorption. 

In Fig. [9l we demonstrate the principle of x-ray pulse 
shaping based on the effect electromagnetically induced 
transparency for x rays. On the left hand side of the gas 
cell, the initial laser and x-ray pulses are shown. The 
laser pulse has the shape of Eq. (|4T|l : the x-ray pulse 
has a Gaussian envelope ([55)1 . The pulses copropagate 
from left to right with the speed of light in vacuum c 
through the gas cell which is filled with argon. Also in the 
cell the propagation speed deviates only minutely from c 
(see the appendix). After passing through the cell — the 
right hand side of the figure — a pulse shape similar to the 
shape of the laser pulse is cut out of the broad x-ray pulse; 
the laser pulse remains unchanged. The fairly elaborate 
pulse shape (j41|) was imprinted approximately on the x- 
ray pulse. In Fig. (TUl we show a close-up of the x-ray 
pulse on the right hand side of Fig. [5] for two different 
intensities of the dressing laser. We see a pronounced 
reduction of the initial peak x-ray flux Jx.o [Eq. 
Due to the nonlinear dependence of the x-ray absorption 
cross section on the laser intensity (Fig.|4]), the resulting 
transmitted pulse varies not only in the flux but also 
the general shape differs somewhat from the initial pulse. 
Overall the peak x-ray flux drops for 10^'^W/cm^ by a 
factor of 10. For neon, we found a stronger transmission 
on resonance which leads to a suppression of only a half 
for the same peak laser intensity p^ . 

Only the flux envelope of the x-ray pulse has been mod- 
ified in Fig. [TOl The time evolution of the phase of the 
x-ray pulses remains random. With the EIT for x rays 
method, the phase cannot be effectively manipulated be- 
cause the magnitude of the polarizability and variation 
of the polarizability due to laser dressing are very small. 
Hence the magnitude and variation of the index of re- 
fraction of the gas is also very small. The x rays are 
basically fully absorbed before the phase can be influ- 
enced noticeably. This finding is in contrast to EIT for 
optical wavelengths where enormous changes of refrac- 
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FIG. 8: (Color online) X-ray polarizability of argon for per- 
pendicular laser and x-ray polarization vectors. The dressing 
laser has an intensity of lO^'^ W/cm^ in the upper panel and 
of 10^'' W/cm^ in the lower panel. 



tion and dispersion of the medium on the resonance can 
be observed and exploited see, e.g., Refs. ll§L fl9l [20l [2ll. 



V. CONCLUSION 

In this paper, we investigated the complex index of re- 
fraction of atoms in the light of an intense optical laser 
in the x-ray region. This is a so-called two-color prob- 
lem. The laser intensity is assumed to be sufficiently low 
such that the atoms are only dressed but not excited or 
ionized. We devised an ab initio theory to compute the 
basic atomic quantities, the dynamic dipole polarizabil- 
ity and the photoabsorption cross section, from which the 
real and imaginary part of the index of refraction, respec- 
tively, are determined. We describe the fundamental in- 
teraction of the atoms with the light in terms of quantum 
electrodynamics. The determination of the atomic prop- 
erties involves resolvents. We use a single-vector Lanc- 
zos algorithm to compute the involved resolvents directly. 
The index of refraction is a concept of classical electrody- 
namics; it follows from the Maxwell equations. The tran- 
sition from the quantum mechanics to classical physics is 
made by expressing the classical macroscopic polarizabil- 
ity in terms of quantum mechanical atomic polarizability 
and atomic absorption cross section. Finally, we consider 
a laser and an x-ray pulse copropagating through a gas 
cell using rate equations. 

Our theory is applied to argon. We study the x-ray po- 
larizability and absorption over a large range of dressing- 
laser intensities for parallel and perpendicular laser and 
x-ray polarization vectors. For parallel polarizations, we 
find electromagnetically induced transparency (EIT) for 
X rays on the Ar Is — > Ap pre-edge resonance. The ab- 
sorption is suppressed by a factor of 3.6 when the laser is 
present. Also for perpendicular polarizations, the cross 
section on resonance shows a noticeable drop by a factor 
of 1.9 for a laser intensity of 10^^ W/cm^. As an applica- 
tion of EIT for X rays, the control of the absorption of the 
X rays allows one to imprint the shape of the laser pulse 
on X rays. We show the transmitted x-ray pulses for two 
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FIG. 9: (Color online) X-ray pulse shaping with electromag- 
netically induced transparency. Propagation of a laser pulse 
[Eq. (|4H) ] (red, high, narrow) and an x-ray pulse [Eq. (|38|l ] 
(black, Gaussian, broad) on the right hand side through a cell 
(green, hatched) filled with argon gas. The pulses copropa- 
gate from left to right with the speed of light. 



different laser intensities which vary, apart from the flux, 
also due to the nonlinear dependence of the cross section 
on the laser intensity. 

Our work opens up numerous possibilities for future 
research. The K edge of argon is at 3205.9 eV [l| which 
implies a wavelength of the x rays of 3.9 A, i.e., our work 
opens up possibilities to shape hard x-ray pulses. Ultra- 
short, femtosecond, hard x-ray pulses come into reach. 
Especially, in conjunction with the emerging x-ray free 
electron lasers, our results may prove useful. The ability 
of pulse shaping x rays opens up possibilities for all-x-ray 
pump-probe experiments and quantum control of inner- 
shell processes. 
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APPENDIX: SECONDARY PHYSICAL 
PROCESSES 

Numerous secondary physical processes occur in course 
of the pulse propagation. The absorption of x rays leads 
to the production of ions and free electrons; in other 
words, a plasma. The ion production is governed by the 
rate equation 



(A.l) 



In the inner part of the gas cell the ion number density 
is constant; its maximum is 1.3 x 10^^ cm^"^. The ions 
which are created by photoionization have a K vacancy. 
Each resonantly absorbed x ray generates a decay cascade 
which leads to the emission of 3.8 electrons (Fig. 1(b) in 
Ref.i61|)- 
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FIG. 10: (Color online) The x-ray pulse shapes that are ob- 
tained in argon gas with the setup in Fig |9] for two regimes of 



laser intensity: 
(lower panel). 



lO^^W/cm^ (upper panel) and lO" W/cm^ 



The phase velocity (53[ of laser and x-ray pulses in the 
medium is 



Wp(w) 



n{ijj) 

and their group velocity is pll . [ssj : 

r 

Vg{uj) 



(A.2) 



(A.3) 



where n{oj) denotes the real part of the refractive index 
in this appendix. We will show that both velocities differ 
only minutely from the corresponding velocities in vac- 
uum for our choice of parameters. 

The index of refraction of a plasma is 



ipi(a;) = 



with the plasma frequency ujp — y/^nne 
velocity for a plasma is 

Wg,pi(w) = cnpi(w) . 



(A.4) 

[H. The group 



(A.5) 



The refractive index of argon gas for laser light is 
obtained from Eq. ((37|) . We use the dynamic polariz- 
ability of argon atoms for wl = 1.55 eV [6^ [6^ which 
is 11.42 flp E/i~^. The real part of the refractive in- 
dex follows with Eq. §7}i to n(1.55eV) = 1.0013 for the 
physical parameters from Sec. IIIII 

The group velocity of x rays in argon gas follows from 
Eq. (|A.3p . We insert the derivative of the index of re- 
fraction letting fi — in Eq. (|37p yielding 



1 + 2^ n# [ai (wx) + ^ (wx)] 



(A.6) 



The polarizability of argon on resonance is taken from the 
curve in Fig. [5] and its numerical derivative. The group 
velocity of the laser in argon gas Vg^i^{uji,) is obtained 
from Eq. (|A.6|) — with "X" replaced by "L" — by inserting 
the slope of the line used for the linear interpolation of 
the dynamic polarizability [6^ . [63| . 
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LO [eV] 


Vp/c - 1 


^g/c- 1 


Vp^pl/c - 1 


^g,pi/c - 1 


1.55 


-1.3 X 10"^ 


-1.3 X 10"^ 


1.4 X 10"* 


-1.4 X 10"* 


3203.42 


9.6 X 10"* 


4.0 X 10~* 


3.4 X 10"^^ 


-3.4 X 10"^^ 



TABLE I: Estimate of the changes of the phase and group velocity due to secondary physical effects. We consider refraction 
and dispersion of the laser and the x rays by the argon-gas medium and the plasma for /l,o = 10^'' W /era? . 



In Table [H we list exemplary estimates of the sec- 
ondary physical effects. We assume a laser intensity 
of lO^^W/cm^. Overall we observe that these effects 
are very small under the conditions of this paper. How- 
ever, for increasing gas densities and/or increasing laser 
intensities, these factors become relevant and have to be 
accounted for. The largest difference to c occurs for the 
phase and group velocities of the laser in the argon gas. 
We observe that the maximum of the phase velocities and 
the group velocities are essentially the same. This is un- 



like EIT for optical wavelengths where the group velocity 
and phase velocity on resonance differ largely, leading to 
slow light [H, [2l|. For lower laser intensities, the plasma 
related velocities, v-p^p\ and Ug^pi increase slightly because 
a larger fraction of the x-ray pulse is absorbed leading to 
a denser plasma. Conversely, these velocities are smaller 
for higher laser intensities as long as the laser intensity 
remains far below the saturation intensity. As soon as 
the laser begins to contribute to the plasma formation, 
I'p.pi and Ug^pi deviate noticeably from c. 
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